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Generating the primordial curvature perturbations in preheating 
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We show that the primordial curvature perturbations may originate not from the quantum fluc- 
tuations of inflaton during inflation but from isocurvature perturbations which are amplified during 
preheating after inflation. We consider a simple preheating model, whose potential is given by 
= 3^0 4 + |s 2( / ,2 X 2 with g 2 /\ — 2, as a possible realization of generating curvature pertur- 
bations during preheating. We make use of the SN formalism which requires only knowledge of the 
homogeneous background solutions in order to evaluate the evolution of curvature perturbations on 
super-horizon scales. We solve the background equations numerically and find that the amplitude 
and the spectral index of curvature perturbations originating from preheating can be tuned to the 
observed values if the isocurvature perturbations at the end of inflation are not suppressed on super- 
horizon scales. We also point out that the tensor to scalar ratio in jXfj) 4 inflation model can be 
significantly lowered, hence letting the |Acjf> 4 model, which is ruled out by the Wilkinson Microwave 
Anisotropy Probe (WMAP) data combined with SDSS data, get back into the observationally al- 
lowed region. 

PACS numbers: 98.80.Bp, 98.80.Cq 



I. INTRODUCTION 

The primordial curvature perturbations which are co- 
herent over the horizon scales at recombination epoch 
are seeds for large scale structure. A simple scenario 
to explain the origin of curvature perturbations is that 
vacuum fluctuations of inflaton driving inflation convert 
into classical perturbations at around the time of hori- 
zon crossing, which are almost adiabatic, almost scale- 
invariant and almost Gaussian. These perturbations are 
constant until the scales of interest re-enter the horizon. 

Recently alternative mechanisms were proposed which 
generate the primordial curvature perturbations after in- 
flation, not from the inflaton fluctuations. In curvaton 
scenario 0, H Q, quantum fluctuations of light scalar 
field during inflation which arc not necessarily related 
to the inflaton, convert into curvature perturbations af- 
ter inflation ends. In inhomogeneous reheating scenario 
Qilll El; spatial fluctuations of the decay rate of the infla- 
ton to radiation generate curvature perturbations. Both 
of these mechanisms are based on the formula 
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which is valid on super-horizon scales. Here ' is a deriva- 
tive with respect to the cosmological time, H = a/a rep- 
resents the Hubble parameter, p is the energy density, P 
is the pressure, 1Z C is the curvature perturbation on co- 
moving hypersurface and <5P n ad = SP — ^-Sp is called the 
non-adiabatic part of the pressure perturbation or isocur- 
vature perturbation. If <5P n ad is non-zero after inflation, 
which is the case of the curvaton and the inhomogeneous 
reheating scenarios, 1Z C on super-horizon scales is gener- 
ated after inflation. 

As another candidate for the mechanism of generating 
curvature perturbations after inflation, preheating model 



has also been considered [1, 01 El El 0- Preheating 0, 
HH ITU El is a process in which the coherent oscillation 
of the inflaton excites the fluctuations of fields coupled to 
the inflaton by parametric resonance. This process was 
investigated in detail E3, EH H3, [S3, E3, H3, 

El Elll El El El 13 without metric perturbations, 
which is valid in dealing processes occurring well within 
the horizon scale. It was then recognized that particle 
production in preheating occurs much more efficiently 
than in the old reheating scenario where inflaton decays 
perturbatively to other particles. 

Naively, fluctuations of such fields coupled to the in- 
flaton seem to be excellent source of isocurvature pertur- 
bations. This means there is a possibility of generating 
curvature perturbations on super-horizon scales during 
preheating epoch. However, whether preheating can af- 
fect the evolution of curvature perturbations on super- 
horizon scales or not [1 0] is not completely understood 
yet. A necessary condition that curvature perturbations 
on super-horizon scales are affected by parametric res- 
onance is that isocurvature perturbations are not sup- 
pressed on super-horizon scales j^lj. Though this condi- 
tion is not satisfied for large class of two-field preheating 
models, there remains some models which satisfy this 
condition. A simple example is a massless self-coupled 
inflaton (f> coupled to another scalar field Xj i- e - with 
potential [H] 



A 



o 2 



(2) 



with g 2 /X = 1-3. 

For the model given in Eq. with g 2 /\ — 2, where 
the k = mode has the largest growth rate(fc is the 
wavenumber), it was suggested in Ref. that large- 
scale curvature perturbations are much amplified and ex- 
ceed the current observational upper bound by using the 
mean field approximation. 



On the other hand, in Ref. [TT|| . by making use of the 
SN formalism [23, 13 HH > it was suggested that the am- 
plification of curvature perturbations during preheating 
is less efficient than is expected in 0, [3~2 . 

In the SN formalism, we can evaluate the evolution 
of curvature perturbations only by the knowledge of the 
homogeneous background solutions and as was empha- 
sized in [ll| , unphysical acausal energy transfer does not 
occur in SN formalism while this energy transfer is not 
inhibited in the mean field approximation used in 

However in Ref. whether the generation of cur- 
vature perturbations on super-horizon scales takes place 
during preheating or not remains an open question be- 
cause the number of solving background equations during 
preheating varying the initial value of the isocurvature 
perturbations at the end of inflation is not enough. 

In this paper, we numerically solve background equa- 
tions for the model given by Eq. Q during preheating 
and show that scale invariant curvature perturbations 
with amplitude observed in the temperature anisotropy 
of Cosmic Microwave Background(CMB) can be gener- 
ated during preheating. 1 

Before we end this section, we define the power spec- 
trum of any field f(x), which we denote as Vf, by 
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Here (■ • • ) is the ensemble average and k is the wave 
number vector. 

The paper is organized as follows. In Sec. II we briefly 
review the SN formalism which we use in a later section in 
order to evaluate curvature perturbations generated dur- 
ing preheating. In Sec. Ill we give the basics of preheat- 
ing and give homogeneous background equations which 
we numerically solve. In Sec. IV we solve the background 
equations numerically and discuss how we can obtain the 
power spectrum of curvature perturbations from the so- 
lutions of the background equations. Then we give the 
resulting power spectrum after preheating. Sec. V is a 
summary. 



II. SN FORMALISM 

In this section, we briefly review the SN formalism [3~2l 

MM 

Let us consider two world lines A and B which repre- 
sent the evolution of two separate homogeneous universes 



1 Generating the primordial curvature perturbations in preheating 
was als o d iscussed in ll j . They considered instant preheating 
model 36] in which case the resulting curvature perturbations 
are proportional to the initial isocurvature perturbations. As we 
will show in Sec. IV, the resulting curvature perturbations for 
the case we consider are very sensitive to the tiny differences of 
the initial isocurvature perturbations. 
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FIG. 1: The schematic picture of the 57V formalism 



separated by the distance d, and take U n i to be a certain 
time soon after the relevant scale d, becomes larger than 
the horizon scale cH , and i nn to be a certain time after 
the complete reheating. 

As is well known, at t > ig n the curvature perturbation 
on comoving hypersurface, 1Z C becomes constant in time 
on super-horizon scales. We define the background e- 
folding number as 



N 



Hdt 



(4) 



According to [33J, |3J, |3Jj , taking the initial (t = t ini ) hy- 
persurface to be flat and the final (t = tfi n ) hypersurface 
to be comoving, there is a simple relation between 1Z C and 
the perturbation of the background e-folding number as 

n cd {t Rn ) = - [N A {tfi n ,t ini )- /V B (ifin,*im)] . ( 5 ) 

where TZ c d represents the comoving curvature perturba- 
tion on super-horizon scale, d, and Na(b) represents the 
e-folding number measured from t; n j to ifi n along the 
world line A(B). Thus, in order to evaluate the final 
curvature perturbation on super-horizon scales, it is only 
necessary to calculate the difference in the e-folding num- 
ber measured from the initial to final hypersurfaces be- 
tween separate FRW universes. It should be also noted 
that we only assume that the metric perturbations on 
super-horizon scales are linear perturbations, while we 
do not assume the linearity of the matter fields in the 
FRW universes here. Hence we can also use this formula 
for the case in which the dynamics of matter fields is 
non-linear, such as in preheating. 

In our calculation, we obtain the e-folding numbers 
of separate homogeneous universes from i; n j to tfi n , given 
the different values of x on initial hypersurface at t = t ml . 
We take the t = t^ n hypersurface to be constant Hubble 
one, that is, Ha — Hb, rather than comoving hypersur- 
face because it is much easier to handle in the numerical 
calculations. It is well known that the constant Hubble 
hypersurface is identical to the comoving hypersurface 
after complete reheating, that is, during the universe is 
characterized by a single parameter, e.g., the energy den- 
sity of the radiation. 
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III. PREHEATING MODEL 
A. Model 

We consider a two-field preheating model whose effec- 
tive potential is given by Eq. (J2J. In the linear stage 
of preheating where fluctuations of both (j) and x have 
negligible effects on the background dynamics, the ho- 
mogeneous part of 4> oscillates under the potential jcj) 4 . 
The explicit form of the oscillation is a -1 <p cn(x; -^), 
where a is the scale factor, 4> Q is the initial amplitude of 
the oscillation of (f>, x = \/\<j) r] (77 is a conformal time) 
and cn is the Jacobi cosine function. In this regime, the 
evolution equation for a given mode k of the x field is 
given by [2f| 

X'i + (k 2 + ^cn 2 (x; X k = 0, (6) 

where = axt, k 2 = k 2 /(\(f> 2 ) and ' denotes the deriva- 
tive with respect to x. It is known that for n(2n — 1) < 
g 2 /X < n(2n + 1) (n is a positive integer), Eq. © has 
a solution which grows exponentially at k = 0, i.e., in 
the long wavelength limit. Interestingly, for g 2 /X = 2, 
the longest wavelength mode has the largest growth rate 
and the power spectrum of x-field at the end of inflation 
is scale- invariant |10| if the interaction between (f> and x 
during inflation is also given by Eq. ©. 

If we allow the non-minimal negative coupling \£,Rx 2 
as an additional interaction to the effective potential 
Eq. © 37], the effective mass squared of x-ficld is 

m 2 «. 9 V+8^ 2 ^) 2 . (7) 

Since the second term decreases faster than the first one 
as <f> decreases, it is possible that the second term domi- 
nates during inflation but becomes negligible during pre- 
heating. In this case, the power spectrum of x a t the end 
of inflation can be red rather than scale- invariant. In 
the following discussions, we do not fix the initial power 
spectrum of x-field and study the resulting power spec- 
trum of the curvature perturbation after preheating for 
various cases of V x . 

At the linear stage of preheating, x-field which corre- 
sponds to the isocurvature perturbation |3l| grows ex- 
ponentially by parametric resonance (</>-field fluctuations 
also grow but its growth rate is much smaller than that 
of x-field. [2^|) Eventually non- linear interactions of fields 
dominate the dynamics. Thus, in order to evaluate cur- 
vature perturbations, we use the SN formalism instead 
of solving the perturbation equations. 

B. Background equations 

In the two-field preheating model in which the effec- 
tive potential is given by Eq. J2J), the homogeneous back- 



ground equations are 

4> + (3H + T)4> + A0 3 + g 2 <\>x 2 = 0, (8) 
X + (3ff + r)x + <?Vx = 0, (9) 
p R + 4Hp R = T(4> 2 + x 2 ), (10) 

H 2 = ^(p S + Pn), (ID 

Ps ^ 2 + ix 2 + ^ 4 + ^V. (12) 

Here pn is the energy density of the radiation. We have 
included the radiation fluid in order to realize the radi- 
ation dominated universe after preheating. T is a decay 
rate from scalar fields to the radiation. We assume that 
both the decay rates of <fi and x to the radiation are the 
same for simplicity. 

IV. NUMERICAL RESULTS 
A. Dependence of 8N on the initial value of \- 

In order to evaluate curvature perturbations generated 
during preheating by using the 5N formalism, we first 
study the dependence of the e-folding number N until 
the Hubble parameter H reaches a given value -fffi n vary- 
ing the initial values of fields, that is, fields values at the 
end of inflation. Perturbation of <f> approximately corre- 
sponds to the adiabatic perturbation in the linear regime 
of preheating. Hence the variation of (f> at the end of in- 
flation does not newly-generate curvature perturbations 
during preheating and we only vary the initial value of x 
in the numerical simulations. 

We show in Fig. |3 plots of e-folding number as a func- 
tion of the initial value of x which we denote as Xini for 
two different values of the sample size of Xini (10 2 and 10 4 
for the left and right panel respectively). The parame- 
ters used in the simulations are g 2 = 2A = 0.04, (j>- m i — 

m p , 4n = xmi = 0, r = io- 4 # ini , H &n = io- 1 r. 

From both plots, we find that the e-folding number N 
is very sensitive to Xini- This feature is also seen in \H\ 
though the sample size is limited to 10 2 and the param- 
eters such as decay rate T arc different from those used 
in our simulations. Hence we expect that this feature 
is characteristic of preheating itself. We have increased 
the sample size to 10 6 (we cannot show the plot in this 
case because the file size becomes too large.) and found 
similar feature to the case of the smaller sample size. 

This sensitive fluctuations of N can be interpreted in 
a following way [TlJ. Neglecting the radiation to which 
both (f> and x decay for simplicity, the total energy density 
is given by ptot = Ps, and the evolution equation of ptot 
is 

^ = -4ptot - H-^ (W,) - 3if<^, (13) 

where i labels <f> and x- In the oscillating phase after infla- 
tion, the second and third terms oscillate. If the period of 
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these oscillations are much smaller than the Hubble time, 
the averaged energy density over the period will decay 
in proportion to e~ 4JV , like the radiation energy density. 
In such a case, the evolution of the universe can be de- 
scribed by only one parameter, such as a radiation energy 
density, and the isocurvature perturbations do not arise. 
Hence additional curvature perturbations are not gener- 
ated during preheating on super- horizon scales. How- 
ever, the sign of (j>i4>i remains the same for a long time 
when the background trajectory is trapped in the poten- 
tial valley along 4> — 0. We plot in Fig.|3]the background 
trajectories of scalar fields for two slightly different initial 
values of x (left panel is Xini = 5.1 x 10~ 9 M p and right 
one is Xini = 5.2 x 10~ 9 M p ) with other parameters fixed 
to the same values. We see that these two trajectories 
soon deviate and each trajectory becomes chaotic [38|. 
In particular, the behavior of 4>i<fii is quite sensitive to 
the tiny difference of Xini- We found that when the back- 
ground trajectory is trapped in the potential valley along 
= 0, which is seen remarkably in the right panel, the 
time interval between two zeros of </>j(/>i becomes so long 
that it is smaller than the Hubble time only by a factor. 
In such a case, the expansion rate of the universe devi- 
ates significantly from that of the radiation dominated 
universe. Because the background trajectory is chaotic, 
the entrapment of the trajectory in the potential valley 
along <j) = occurs almost randomly for Xini- Hence the 
dependence of N on Xini becomes almost random. 

Before we go to the discussion of how to obtain the 
power spectrum of curvature perturbations from the data 
which gives the correspondence between Xini and the e- 
folding number AT, we show in Fig. ^ the dependence 
of the mean and the variance of N on the sample size, 
which we denote as p. From this result, we find that 
the dependence of both (N(x)) x and (SN(x) 2 ) x on p be- 
comes weaker as p becomes larger. Here 5N(x) is defined 



by 

SN( X ) = N(x) - (N( X )) X , (14) 

and (■■■) x is the average over the sampling points of 
Xini- In particular, the fluctuations of (5N(x))x f° r 
different values of p are within the statistical ones ~ 
v/ (SN(x) 2 ) x /p, which suggests that both the mean and 
the variance converge in proportion to 1/y/p for p — > oo 
if we fix the range of x in which we take average. Hence 
the smoothed 6N(x) over Xini defined as 

SN(A X ;x) =J2 S * w x(&r,Xj-x) SN( Xj ), (15) 

Xj 

also converges if the large number of sampling points are 
contained in the range (x — A%,x + Ax). Here Xj is 
the sampling point (integer j labels the sampling point) 
and 5x = Xj+i ~ Xj- W x (Ax; x) is the window function 
which rapidly drops for |x| ^ Ax and its normalization 
is defined by 

[d x W x (A X ;x-x') = l- (16) 

In Fig. |21 we see that the standard deviation of 6N(x) 
which is about 2 x 10~ 4 is larger than the COBE normal- 
ization 331 . But this does not mean that curvature per- 
turbations on large scales exceed COBE normalization. 
What we actually observe is the perturbation of e-folding 
number smoothed over many different spatial points in 
the real space, or equivalently, average of SN(x) in Xini 
space with some weight factor. In the next subsection, 
we will show the relation between SN(x) and V^^^k). 
Here we denote the perturbation of the e-folding num- 
ber at a point x in the real space as SN(x) rather than 
6N(x) in order to avoid confusion because we use 5N as 
the perturbation in x-field space. 
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FIG. 3: Plots of the background trajectories of and x during preheating. In left and right panel, Xini is 5.1 x 10~ M p and 
5.2 x 10 _9 M P with other parameters fixed to the same values. We find that two trajectories soon deviate and become chaotic. 
In left panel, the background trajectory is trapped in the valley <f> — for a long time. 
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FIG. 4: Plots showing the dependence of the mean (left panel) and the variance (right panel) of N on the sample size p. We 
see that both the mean and the variance converge for large sample size. 



As we mentioned above, SN(x) is very sensitive to \ 
and the dependence of 5N(x) on X seems almost random. 
If this dependence is completely random, then SN(x) is 
also random regardless of the shape of the power spec- 
trum of x & t the end of inflation. Hence, typical magni- 
tude of 5N(x) smoothed over the scale R which is beyond 
the horizon scale at the time of preheating is given by the 
statistical fluctuations as 

~y/Wttf) x (£) ^ (i7) 

where tr is the horizon scale at the time of preheating. 
Hence the smoothed SN(x) is significantly suppressed 
and the generation of curvature perturbations on super- 
horizon scales does not occur. On the other hand, if 
5N(x) is n °t completely randomized, then the amplitude 
of the smoothed SN(x) would be larger than that of the 



statistical fluctuations, in which case it may be possible 
that curvature perturbations can be larger or comparable 
to the COBE normalization. 

In order to see whether SN(x) is randomized or not, 
we divided the region of Xini which contains one million 
sampling points into ten bins and took the average of 
SN(x) for each bin with the same parameters as in Fig. El 
Figure|S]is the result. We see that the fluctuations exceed 
10~ 5 by order of magnitude which is much larger than 
what is expected if SN(x) is random. Hence SN(x) is 
not completely randomized. 

B. Construction of the Power Spectrum of 8N 

Generally, the power spectrum of the completely ran- 
domized field is proportional to fc 3 , that is, extremely 
blue |40j. In the previous subsection, we have found that 
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FIG. 5: Plot showing the e-folding number TV averaged in 
each bin in which 10 sampling points are contained. We 
find that the fluctuation of the averaged e-folding number 
exceeds the statistical fluctuation. Hence TV is not completely 
randomized. 



SN(x) is not completely randomized, so we expect that 
Vgftik) differs from the extremely blue power spectrum 
for sufficiently large scales if the isocurvature perturba- 
tions at the end of inflation, which correspond to x-field 
perturbations in this model, are not suppressed on large 
scales. To see this, we first define cr(R) as the mean of 
the standard deviation of x i n the region of size i?, 

-i 

a 2 (R) = (( X (x)-x(R;x)) 2 ) R ^ jT* fr x (k), (18) 
where x(fy x) 1S defined by 

Here W(R; x' — x) is a window function which rapidly 
drops for \x' — x\ > R and its normalization is given by 

J*.W(KS-t)-X. (20, 

We also define ^(R) as the standard deviation of \ 
smoothed over the scale R, 

r°° dk ? 

^(R) ee (x 2 (R;x)) =J liT \W(R- l ;k)\ V x (k), 

(21) 

where W(R 1 ;fc) is a Fourier window function and we 
have introduced a box of size L(3> R), within which the 
stochastic properties are to be defined. 

A key quantity for calculating the power spectrum of 
curvature perturbations on super-horizon scales is SN(x) 
smoothed over the scale R> r#, 

SN(R;x) = r 3 H J2 W ( R ^i ~ ^)( N ( £ i) ~ (%), ( 22 ) 

Xi 



where we have discretized the space into lattice whose 
grid spacing is and have replaced the integral over the 
space with the sum over the lattice points xl in the box. 
(■ ■ ■) l denotes the expectation values inside the box. In 
Eq. 122fl . we implicitly assume that the spatial volume of 
size R contains sufficient large number of lattice points 
so that the contamination of the random component is 
significantly reduced and SN(R;x) does not depend on 
the size of m- This assumption is valid as long as we 
consider sufficiently large scales (but R <C L is imposed 
through). 

In order to relate 5N(R; x) to 6N(x) averaged over the 
sampling points of Xini, we rewrite Eq. I|22|l as the sum 
over the sampling points on Xini space, 

5N(R; x) = ( ! |) 3 E N <Xi) S i Xi) - (n(R, x; Xj ))l) , 

X] 

(23) 

where we have used the top hat window function. 
n(R, x; Xj) is a number of spatial points in the region of 
size R with x at its center in which Xini takes the value 
between Xj an d Xj+i- From the definition of n{R 1 x; Xi)i 
we have the following relation, 

Y,n{R,x; Xj )= (— ) 3 - (24) 

From this relation, we can add the Xj -independent quan- 
tity to N(xj) in Eq. I|23() . So, as a matter of convenience, 
we can rewrite Eq. I|23|) as 

8N(R; x) = (^) 3 J2 W(Xi) % Xj) ~ (n(R, 2\ Xj))l) . 

Xj 

(25) 

Eq. 1)2 5[1. or equivalently Eq. (J2SJ, is not strictly correct 
because from Eq. (1221) to (|25|l we have replaced x{xi) with 
Xj which satisfies Xj < x(xi) < Xj+i an d this will induce 
an error in the corresponding replacement of 5N{x) of 
order */ (5N 2 ) X for each lattice point. The sum of this 
error of SN(x) over the lattice points would roughly cor- 
respond to the variation of (SN(x))x wnen each sampling 
point of Xini is shifted by 0(6x)- We have verified that 
the variation of (SN(x))x under the shift of Xini 

by 0(S X ) 

lies within (SN{x) 2 ) x / ' y/p. Hence Eq. ggj is correct 
to good accuracy as long as R^> rjj. 

If 5x is small enough, n(R,x;x) will be proportional 
to 5x- Hence we introduce f(R, x; x) as 

n(R, x- x ) = (J£) f(R, x; X )S X - (26) 

Then f(R, x; x) satisfies the following normalization con- 
dition, 

J d X f(R,x;x) = l- (27) 
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Using Eq. (j^SJ, we have 
(SN 2 (R;x)) L = ^^(<5 X ) 2 SN( Xj )SN( Xk ) 

Xj Xk 

x (f(R,x;xk) - (/(i?,x;x fc )) L )) L . 

(28) 

We have to specify the function form of f(R, x; x) which 
depends on the model in order to proceed further from 
Eq. J2HI- We assume that f(R,x;x) takes the form as, 



Substituting Eqs. ^ and (E3 into Eq. pfll . we have 
(8N 2 (R;x)) L = EE (^) 2 5N (Xj)SN(Xk) 



Xj Xk 



i! (J?0(cr^(J?,) + 2n i! (H)) 



27rcr(i?) y/a 2 (R) + 2fi 2 (R) 



2 i 2 



2 7 r(/i 2 ( J R) + a 2 (iZ)) 



(34) 



f(R,x; X ) 



1 



27rcr(i?) 



exp 



2(T 2 {R) 



(29) 

This assumption is valid if x-field obeys the Gaussian 
distribution, which is the case if x-held is in the vacuum 
state at the end of inflation. Under this assumption, we 
have 



Eq. I|34(l relates the power spectrum of curvature pertur- 
bations to the average of 5N(xj)SN(xk) over Xj an d Xk 
with some weight factor. 

If the power spectrum of x-field at the end of inflation 
is red, we can reduce Eq. 134|) to a simpler form. In 
this case, we have /i(-R) ^S> a(R). Then Eq. (|34|1 can be 
approximated as 

{8N 2 (R; x)) L ^^ (Sxf SN( Xj )5N( X k) 



Xj Xk 



{f{R ^ ))L = V^iR~) 

{f(R,x; X j)f(R,x;Xk))L = 



(e -^yfe-x(^» 2 )Lj (30 ) 



2V2ira(R)fi(R) 



2ira 2 (R) 



( e -^Hy^-x(fl;S)) 2 



27T/X 2 (i?) 



e 2 ^ 



■^(x ] -Xk) 2 -j^(x 2 J +xl) 



xe 2CT ^ fl > 



^(xk-x(R;S)) 2 



(31) 



d X 



-oo " (2irf /2 a 2 {R)Li{R) 



]T$> X 2 SN( Xj )SN( Xk 



Xj Xk 



We could replace ensemble average in Eqs. (|3lHl and (fTTTf) 
with the spatial average in the box of size L. Because x 
appears in Eqs. l|3"U|) and l(3"TJl only through x(i?;af), the 
ensemble average in these equations can be also inter- 
preted as the ensemble average over x(R'i ^0 with some 
weight factor. If %-field is the Gaussian random field, 
then x(R'i x) is also Gaussian random field, in which 
case the ensemble average of Eqs. <|3(J[) and <|31[) becomes 
Gaussian integral over x(R'i Then we have 



xe 



E** 



L Xj 



2tth{R) 



dxSN 2 (a(R);x) 



G(a(R),n(R)). 



1 



2tt(i(R) 



exp 



2,i 2 (i?) 



(35) 



1 



2irfi(R) 
a(R) 



dx(R;x)e 2 ^<«> e 2 



M##„-^(x 3 -x(^)) 2 



y/a*(R)+IJ?(R) 



(32) 



and 



/g 2<^(fl) 1 
1 



2 7 T J u(i?) J- 

xe -^(«)fe-x(^)) 

a(i?) 



e -i^)te-x(^)) 2 x 



<2%(-R; f) e 2 ^ 2 (-R) 



e -^yte-x(^)) 2 



M 2 (A) 



/ 2 2x2 

r(x,-Xfc) 



V^(i2) + 2^(R) 



(33) 



Hence (8N 2 (R;x))l is equal to the variance of 
5N(a(R);x) m the range \x\ ^5 A*(-R)- We expect that 
V s fj{k) shows deviation from the extremely blue power 
spectrum for scales k < i?" 1 where R c is the critical 
scale at which the fluctuation due to the randomness 
~ ( ~gy (<^(x) 2 )x S oes down to the value given by 

On the other hand, we have 

^ 25/*" = 

(36) 

From this simple relation between (5N 2 (R; x))l and the 
power spectrum of SN(x), we can make a connection be- 
tween "P^ and 5N 2 (a(R); x) 

given by Eq. I|35l) . 

Furthermore, in order to have an analytic relation be- 
tween Vgfj and 8N 2 (a(R);x)i we assume that V x (k) is 
given by 



n x > 0. (37) 
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Then fi(R) and cr{R) become 

fi 2 (R)~ — V (k L) nx , a 2 {R)^—V {k R) n \ (38) 

where L is the size of the box as mentioned before. We see 
that does not depend on R and a(R) is an increasing 
function with respect to R. Using Eq. 1)38(1 . we have 

From Eqs. (|39|) and (|36[) . we obtain the relation between 
V sf4 and G(a(R),n(R)), 

-PstiiR- 1 ) = -Y aiR) d^R) G(a(jR) ' (40) 

Hence we can obtain the power spectrum of curvature 
perturbations at any scale if we know the dependence of 
a{R) on G(a(R),fi(R)). 

The top left panel in Fig. plots G(a(R), n(R)) as a 
function of cr(R) for fixed n(R). In this plot, we have 
approximated G(a(R), £t(i?)) as 



G(a(R),fi(R)) 



MR) 



d X 6N 2 (a(R); X ), (41) 



for simplicity. We assume that this replacement for cal- 
culating G(a(R), fi(R)) does not cause significant error 
on the evaluation of the amplitude of curvature pertur- 
bations. n(R) is chosen to be 5 x 10 -5 M p . Solid line 
corresponds to the case of 10 6 sampling points of Xini 
within the region < Xini < l^(R) and dotted one corre- 
sponds to the case of 10 5 sampling points. 

We first see that slopes of both two lines become steep 
at small a{R). This steep rise is due to the contamination 
of the random component. Since the sampling number in 
the width of fJ-(R) is limited in our simulations (at most 
10 6 ), the sampling number in the width of cr(R) becomes 
at most ~ 10 6 x ^^y- Hence if we decrease cr(R) with 
)J,(R) fixed, the sampling number within the width of 
(j(R) also decreases. Below some critical <r(i?), the sam- 
pling number is so small that 5N(a(R); x) is still domi- 
nated by the random component. In Fig. EI we actually 
observe that if we increase the number of sampling points 
the amplitude of G(cr(R), fi(R)) at small <r(R), where the 
slope of the line is steep, decreases and the difference be- 
tween two lines becomes negligible for large a(R). Hence 
this steep rise at small <j(R) is meaningless in a sense that 
we cannot obtain the correct amplitude of curvature per- 
turbations at the corresponding scale R. We have noth- 
ing left except to increase the sampling number in order 
to reduce the contribution from the random component 
and obtain the correct G(a(R), fJ>(R}) at small <r(R). On 
the other hand, at large cr(R), we also see that the slopes 
become steep. However, this feature cannot be trusted 
because the approximation we made in Eq. (|35|l is no 
longer valid for a(R)/fi(R) ~ Oil). In such a case, we 



have to go back to Eq. I|34|l . which requires a heavy nu- 
merical calculation. 

Thus, in the following discussion, when we discuss the 
property of the power spectrum, we consider only for the 
scales corresponding to a(R)/fi(R) ~ e -5 — e -3 . 

The remaining three panels (top right, bottom left and 
bottom right) plot —a(R) a 3 R \ G(a(R), (J,(R)) which is 
proportional to V s fr for three different values of /i(i?) 
with other parameters such as g 2 /A fixed to the same as 
in Fig. [21 jtx(-R) in the top right panel is the same as in 
the top left panel and the solid and the dotted lines in 
the right panel just correspond to the same type of lines 
in the left panel. fi(R) in the bottom left and right panels 
are 10 _8 M p and 10~ 2 M p respectively. We find that the 
amplitude and shape of the plot depend on [i(R). In the 
top right panel, — cr(R) d 3 R \ G(a(R), fJ,{R)) remains al- 
most constant except for the region of small cr(R) where 
the slope becomes steep due to the small size of the sam- 
pling number and the region where a(R)//j.(R) — 0(1). 
On the other hand, —cr(R) d 3 R \ G(a(R), fj,(R)) isamono- 
tonically decreasing function in the bottom left panel and 
monotonically increasing function in the bottom right 
panel. 

Later, we will give a qualitative explanation why these 
plots depend on /i(i?). 

Now let us consider what type of the power spectrum 
we can obtain from the numerical results above. Since 
(j(R) is an increasing function with respect to R, going 
to larger scales means increasing a(R). From Eq. 138fl . 
the horizontal axis of all the panels in Figgis 



In 



°(R) 

KR) 



n v , R 
= -± In — . 
2 L 



(42) 



Hence if we increase the scale of interest R by a fac- 
tor of 10, then the corresponding point on the horizon- 
tal axis shifts by ~ n x . Since cosmologically interest- 
ing scales range from 10 4 Mpc to IMpc, the correspond- 
ing range in the horizontal axis, which we denote as 
A In is ~ 4n x . From Eq. gUJ), the amplitude of V s ft 
is determined by multiplying —o{R) dtJ 9 ^ G{a{ R), n(R)) 



by 



Therefore by choosing suitable n x , we can set 



8N 



to P-r, w 2 x 10 _J corresponding to the COBE 
normalization. For example, if fi(R) — 5 x 10 _5 M p , 
which corresponds to the case of the top right panel, 
—a(R) d J^ R j G(<j(R), fJi(R)) is about in the range 

-5 < \na{R)/ii(R) < -3. Hence if n x w 0.4, then 



V & f] becomes 2 x 10 



9 , which corresponds to the COBE 

1.6. 



normalization. In this case, we have also A In 
If we assume that L has a value so that the horizon scale 



10 Mpc corresponds to In 



-3, then Aln4#T 



lies in the range —5 < \na(R)/ fi(R) < —3, in which case 
V s fj becomes scale invariant over the cosmological scales. 
In the same way as above, for the bottom right panel, if 



n x w 4 x 10 2 , then V S] ^ gives the COBE normalization. 
In this case, A In 



0.2. Because V s <r changes little 
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FIG. 6: In the top left panel, G(a(R), (i(R)) is plotted as a function of a(R). Each solid and dotted line corresponds to the 
case of 10 5 and 10 6 sampling points respectively. The parameters are the same as used in Fig. [5] The remaining three panels 
plot —a(R) G(a(R), fi(R)) as a function of a(R) with fJ,(R) fixed, hence showing the scale dependence of V s ^. /J-(R) in 
the top right, bottom left and bottom right is 5 x 10 _5 M P , 10 -8 M P and 10 -2 M P , respectively. 



in this interval, the spectrum becomes almost scale in- 
variant over the cosmological scales. For the bottom left 
panel, if we choose n x so that P s ^ gives the COBE nor- 
malization, Vsfj becomes too blue over the cosmological 
scales and hence this case is excluded. 

From Fig. we see that there is a tendency that the 
spectrum of the curvature perturbation becomes blue as 
n{R) decreases. This can be explained if we assume that 
5N(x) becomes more random as (J-(R) becomes smaller, 
because the power spectrum of the completely random- 
ized field is extremely blue as mentioned before |40| . In 
our numerical calculation, we devided the range (0, fJ-(R)) 
in Xini space into ten bins and took the average of 6N(x) 
for each bin (cf. see the Fig. [SJ. Then we verified that 
the deviation of fluctuation of averaged SN(x) in each bin 
from the statistical fluctuation becomes smaller as n(R) 
becomes smaller. This result shows that the assumption 
above is valid. 

It seems unnatural that both the amplitude and the 
spectral index of Vgft depend on the size of the box L 
within which stochastic properties are defined. Introduc- 
ing the infrared cutoff L is necessary, because without 



such a cutoff the expectation value of \ 2 diverges. We 
have assumed throughout the paper that the mean value 
of x vanishes in the box. Let us now prepare a box whose 
size is much smaller than L but is still larger than the 
scales of interest inside the box of size L. Then the aver- 
age of x inside the small box does not vanish in general. 
As shown in [ill |42| , statistical quantity such as power 
spectrum defined in the small box is not the same as the 
one defined in the large box of size L. It is only when 
one averages over the small boxes inside the larger box 
of size L that the statistical quantities become the same 
as the one defined in the large box. Hence L, the size of 
the box within which the average of x vanishes, appears 
as a parameter which in this case controls the shape of 

Testing non-Gaussianity of the primordial fluctuations 
has become one of the most powerful probe that distin- 
guishes various inflationary and post-inflationary mod- 
els Recently the level of non-Gaussianity generated 
du ring preheating was investigated in |44j|. The authors 
in |44| consider the same preheating model as in this 
paper. However in |44|. the primordial curvature per- 
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turbations are assumed to be generated during inflation 
while second order fluctuations which are responsible for 
the large non-Gaussianity are amplified during preheat- 
ing. Though it is shown in [44| that massless preheating 
model considered in this paper is already ruled out from 
observations due to large non-Gaussianity, we have to 
reanalyze the problem because the mechanism of gener- 
ating the primordial curvature perturbations in 44] is 
completely different from that in our paper. We are now 
investigating the level of non-Gaussianity generated dur- 
ing preheating by using the 6N formalism. 

Interestingly, the dynamics of preheating is completely 
determined by g 2 /X, not by A and g independently. All 
numerical results shown above is obtained for any value of 
A if we assume g 2 /X = 2. Since A determines the energy 
scale of inflation, A fixes the amplitudes of the primordial 
gravitational waves produced during inflation. On the 
other hand, as we have shown above, if the isocurvature 
perturbations on super-horizon scales are amplified at the 
end of inflation with suitable negative spectral index and 
infrared cutoff, the primordial curvature perturbations 
with almost scale invariant and ~ 10~ 5 amplitude can 
be generated during preheating. Hence we can lower the 
tensor to scalar ratio by lowering A, that is, the energy 
scale of inflation. In this case, A</> 4 inflation model, which 
is ruled out from the WMAP three year results combined 

with SDSS IHElEllilllll, comcs into the observa- 
tionally allowed region. 2 Lowering the tensor to scalar 
ratio in the chaotic inflationary scenario was also con- 
sidered in [§l|(see also H2), though the mechanism is 
completely different from the one in this paper. 

For completeness, we consider the case in which the 
power spectrum of x at the end of inflation is blue, that is, 
fJ,(R) -C a(R). In this case, Eq. (|34[) can be approximated 
as 

(SN 2 (R;x)) ^Y,T, S ^ SN ^ SN ^ 

Xj Xk 

1 1 /2,2\r M 2 (^)/ \2 1 

2wa 2 (R) L J' 

(43) 

Because [xj ~ Xk) 2 is much smaller than unity, 

right hand side of Eq. (|43|l is suppressed. Hence the gen- 
eration of curvature perturbations on large scales does 
not occur in this case. In the scale-invariant case, that 
is, /i(-R) — cr(i?), we have no appropriate approximation 

2 Although, see also loll . 



in order to reduce Eq. i|34[) to a simpler form. So, in this 
case, we have to calculate Eq. (|34f) directly. 

V. CONCLUSION 

The main purpose of this paper is to examine whether 
the generation of the curvature perturbations occurs or 
not during preheating and, if it occurs, whether they 
could be seeds for large-scale structure in the universe. 

To do so, we applied SN formalism to two-field pre- 
heating model whose potential is given by j(f> 4 + \<^X' ''■ 
Then we have analyzed the dependence of the e-folding 
number during preheating on the isocurvature perturba- 
tions at the end of inflation which corresponds to pertur- 
bations of x-field. We verified that the e-folding number 
during preheating is very sensitive to the initial value of 
X due to the chaotic nature of the background trajec tory 
in field space, which has been already observed in |LJ| . 
Though the e-folding number seems almost random as 
a function of Xinii we found that the fluctuation of the 
smoothed e-folding number over a large number of initial 
values of x exceeds the statistical fluctuation. Hence the 
e-folding number is not completely randomized. 

This deviation of the smoothed e-folding number from 
the statistical fluctuation provides the generation of cur- 
vature perturbations on large scales. We numerically 
found that the amplitude and the spectral index of the 
curvature perturbations originating from preheating at 
current cosmological scales can be fixed to the observed 
values by tuning the spectral index of x-field and the 
infrared cutoff. An interesting point is that this conse- 
quence is derived only from the requirement g 2 /A = 2. 
Hence by lowering A which determines the energy scale 
of inflation, we can lower the tensor to scalar ratio to the 
observationally allowed value. 
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